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CONTACT AND CRACK PROBLEMS 
FOR AN ELASTIC WEDGE* 

by 

F. Erdogan and G.D. Gupta 
Lehigh Univ ersity, Bethlehem, Pa. 


Abstract. In this paper the contact and the crack problems for 
an elastic wedge of arbitrary angle are considered. The problem 
is reduced to a singular integral equation which, in the general 
case, may have a generalized Cauchy kernel. The singularities 
under the stamp as well as at the wedge apex are studied and the 
relevant stress intensity factors are defined. The problem is 
solved for various wedge geometries and loading conditions. The 
results may be applicable to certain foundation problems and to 
crack problems in symmetrically loaded wedges in which cracks 
may initiate from the apex. 


1. INTRODUCTION 

In conventional contact problems it is generally assumed 
that- the substrate consists of an elastic half, space or a layered 
medium (see [1] for a thorough discussion). With the application 

in foundation engineering in mind, in these problems the main in- 

\ 

terest has been mostly in the evaluation of the contact pressure. 
If the external load is not applied symmetrically, it is also 
possible to evaluate the ,l tilt angle" of the rigid stamp simu- 
lating the structure [1,2]. In practice one may also encounter 
a certain group of foundation problems in which because of the 
nonuniform stiffness of the substrate the structure may again 
tilt even if the loads are symmetr i cal ly distributed. Among these 
nonsymmetric foundation problems perhaps the simplest one is the 
frictionless contact problem for a plane elastic wedge-shaped sub- 
strate (Figure la). Here the main questions of practical interest 
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are the contact stress distribution and the value of the possible 
tilt angle for a given applied load. In particular one may be in- 
terested in finding the location of the applied load for which the 
stamp (or the structure) would remain vertical, or the angle of ro- 
tation of the stamp if the load is symmetrically applied. 

It may be noted that there is a group of crack problems for 
plane elastic wedges (Figure lb) the formulation of which is iden- 
tical to that of the contact problems (Figure la). Here the cracked 
wedge is loaded perpendicular to and away from the net section 
0=0, b<r<c and with the application to fracture problems in mind 
the main practical interest in the problem is in finding the stress 
intensity factors and the net section stress. In contact as well 
as in crack problems the case of b=0 (Figure 1) is of special in- 
terest. Analytically the problem requires special attention due 
to the fact that for b>0 the formulation leads to a singular in- 
tegral equation with a simple Cauchy singularity whereas for b=0 
the kernel is of the generalized Cauchy type. In the crack prob- 
lem (Figure lb) if the stress intensity factor at b is negative 
the crack located on 0<r<b, 9=0 will close and has to be ignored 
in the solution. In the rectangular stamp problem shown in Fig- 
ure la if the strength of the stress singularity at b is positive 
(i.e.» if the contact stress becomes "tensile") the problem be- 
comes one of "receding contact" [33 with the contact area being 
also an unknown (see [4] for a similar phenomenon in the axisym- 
metric double contact problems). 

In this paper the frictionless contact problem for a plane 
elastic wedge is formulated for an arbitrary stamp profile. As 
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examples the problems for the rectangular (Figure la) and the 
semicircular (insert in Figure 4a) stamps are considered. The 
solution of the crack problem (Figure lb) is obtained by simply 
rei nterpreti ng the results found for the rectangular stamp. 

2 . FORMULATION OF THE PROBLEM . 

Consider a plane elastic wedge of arbitrary angle 0 . Let 
the wedge be subjected to an external load P applied through a 
frictionless rigid stamp of known profile (Figures la). The solu- 
tion of the problem may be obtained by considering either the re- 
lated biharmonic equation or the Navier's equations under the 


following boundary conditions: 

tfge T °» a r e = ( 0 = e 0 » 0 <r<°°) , (l.a,b) 

<? re = 0 , ( 0=0 , 0 <r<“) , ( 2 ) 

= 0, (0 = 0, 0<r<b , c<r<®) , (3,.a) 

|y u 0 = g ( r ) , (9 = 0, b<r<c ) (3.b) 

/ b C c 7 e0 (r ,0)dr = - P, (4) 


where g(r) is a known function and P is. the resultant load per unit 
thickness. Let a and v refer to the following complex stress and 
displacement combinations: 

°(r,e) = a r0 + ia ee ’ 

v(r,0) = |^(u r + i U e ) 5 (0iQl9 o > 0<r<~) . (5.a,b) 
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Using the Mel 1 i n transform one may easily obtain (see, e,g.[5]) 


M[r 2 o] = 2i (s + 1 ) [Ase 1 s9 + B(s+1 )e n ( s+2 ) 0 . Be” 1 ^ s+Z * S ] 

M [r 2 v] = -’2±i[Ase ise + B( s + 1 ) e 1 (s+2) 6 + KBe' 1 (s+2)6 ] , (6.a,b) 

where A ( s ) and B(s) are unknown complex functions and the Mellin 
transform and its inversion are defined by 


M[F) = / “F(r)r s-1 dr, 

r “ n d n F s-1 Hl . _ , ,\n r(s+n) Ilrc , 

/o dr" dr - (_1) ■ rfl J W[F3 ’ 

f (f) = M[F]r- s ds. (7.a-c) 

provided the integral in (7a) is convergent and in (7c) the strip 
of regularity containing the constant y is selected in such a way 
that (the physically imposed) regularity conditions at r-0 and at 
r =co are satisfied. In (6) p and k are the elastic constants of 
the wedge (p is the shear modulus, k=3-4v for plane strain and 
k=(3-v)/(1+v) for plane stress, v being the Poisson's ratio). 


To formulate the problem it is convenient to obtain first 
various Green's functions by replacing the mixed boundary condi- 
tions (3a) and (3b) by the following concentrated force: 

c 00 (r,O) = f 6 ( r-t ) . (8) 


Substituting from (1) 
ft s+1 

sA(s) = Hfs+T) ' 


, (2) and (8) into (6) we find 
sb-j(s) - i(s+2)b 2 (s). 
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B ( s ) = b 1 (s) + ib 2 (s) , 

fJ .s+1 (s+1 ) ( 1 -cos26 )+l -cos2(s+l )e 
b l< s ) ■ 4( s+T) DTD ' 

ft s+] (s+1 )sin2e Q +sin2(s+l >0 O 
b 2 (s) = 4TI+T) Dlil ’ 

0(s) = (s+1 ) 2 (1-cos26 0 ) - [l-cos2(s+l)9 0 ] . (9.a-e) 

If the density function f=f(t) is given or determined, substituting 
from (9) into (6) and inverting, one may obtain the complete solu- 
tion of the problem. In particular, from (5b), (6b) and (9) it 
follows that 

M[r 2 (s + l)b 2 (s). (10) 

Similarly 


M [r 


2 a] = 2i (s + 1 ) [fe ls0 + 2i ( b 1 +i b 2 ) ( s+1 ) e 1 ( S+1 } 0 si he 


+ 2i (b^ -i b^) e“ ^ 9 si n{ s+1 )6] . 


Instead of (8) if we now assume that 

f(r), ( b <r<c ) , 


^ e6 (r,0) = { 


0 , ( 0<r<b , c<r< TO ) , 


OD 


( 12 ) 


from (3b), (10), (9d), (7c) and (12) we obtain the following in- 
tegral equation to determine the unknown function f(r): 

, y+i oo s + 1 

rg(r) = / b c f(t)dt ^ / Y . ico ( ? ) * 


(s+1 )si n26 +sin2(s+l )e 


- ds , (b<r<c ) . 


dTTJ 

In (13) the strip of regularity containing y is determined 
from the behavior of stresses or displacement derivatives as r^O 


( 13 ) 
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and r~>°° , namely that they vanish at r s » and at most may have an 
integrable singularity at r=0. Let 


D(s k ) = 0, (k = +1 , +2,. . . ) , 

Re(s j ._ 1 ) < Re(s - j )<y < Re(Sj )<Re(Sj + i ) , (j=l,2,...). 04) 


Thus, one way to evaluate the kernel in (13) may be to sum the 
residues at s^, (k=l,2,...) for r>t and at s _ k > (k=l,2,..) for 
r<t. Noting that (13) gives Su^/Sr along the infinite line 0=0, 
0<r<«, and the dominant terms for 3Ug/9r are of the order 


■r" (s_i +2 ) ( r< t) , 


~ ( s , t +2 ) 


, ( r>t ) , 


(15) 


the regularity conditions, i.e., 

^ ~r'“, (u<l )for r+O; yp- ~r"“o , (^>1 )f or r-«° (16) 

require that Re(s_^)<-1 and Re( s + i)>-l. If one also notes that -1 
is a root of D{s), it is then clear that the strip of regularity 
of the inversion integrals for stresses and the displacement deriv- 
atives will be 


Re ( s _ i ) < Re ( s ) = y<-l = Re(s + ^) (17) 

where s ^ is the first root of D(s) to the left of the line Re(s)=-1. 

In this problem it is not possible to obtain the roots in 
the closed form. Hence, the resulting infinite series giving the 
kernel cannot be properly studied for the nature and separation of 
possible singular parts. To investigate the singular behavior of 
the kernel it is more convenient to express the inner integral in 
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(13) in terms of a real Integral by letting y=-l indenting the 
contour to the left and defining 

s'+l = iy, log(t/r) = p. (18) 


The integral equation (13) may then be expressed as 
4 ^ - 20_+sin26 


7^7 rg(r) = / f) c f(t)dt[ 


2(26 - 1 +cos 2 0 ) 

O 0 


{ y s i n 2 0 +sinh20V)sinpy 

0 0 dy] , 


. 1 ; “ 

n 0 cosh29^y-l -y^(1 -cos28 o ) 


( b<r<c ) . 


0 9) 


As p-*0 or t/r the infinite integral in (19) becomes divergent. 

Since the integrand is bounded and continuous everywhere in 0<y<~, 
the divergence will be due to the behavior of the integrand at in- 
finity and the divergent part of the integral may easily be sepa- 
rated by considering the asymptotic behavior of the integrand as 
y->«. Thus, adding and subtracting the asymptotic part of the in- 
tegrand to and from the integrand and evaluating the related integ- 
ral, (19) becomes 


4u 

1 + K 


rg(r) = - 1 / fa c f(t)dt[l- - 


tt ( 2 6 +sin29 ) 

0 0 

2(26 2 -l+cos 20 ) 

O 0 


ysin 20 +sinh 20 y 
+ / 0 0 


cosh 2 e y- l -y^ ( 1 -cos26 ) 
o o 


- 1 )si npydy 3 , 


( b<r<c ) . ( 20 ) 


It may easily be shown that for b+0 (20) is an ordinary singular 
integral equation with a simple Cauchy singularity. This may be 
seen by observing that 
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1 - _J ' r1+ f H) n ft . 

rlog(t/r) ,t _ 1) l •, n+1 r 

\ r l 

- ^ [1+0(7 - 1)] - 

In the special case of 0 q -tt 5 using the relation 

/ “(coshay-l )sinpydy = l 2e" 2a ^ n+1 )y sinpydy 
0 n=0 


( 21 ) 


= Z 


Is. 


2 = - l + h COth ia 


( 22 ) 


m=l p^+(2ma) 
the inner integral in (20) becomes 


C ( cosh2u y -T - sinpydy = 7 COth 2 • P • 


(23) 


cosh 27ry^T " 1 ' 2 2 p 

Noting that p=log(t/r) and substituting from (23), it may easily 
be seen that (20) reduces to the integral equation for the elastic 
half plane given by 

^g(r) - - K C ^7 idt ’ (b<r<c) * 


(24) 


which has a closed form solution for any Holder-continuous g(r) [6] 

In another special case where 0 o = 2tt, (i.e., the plane with a 
semi-infinite crack loaded on one flank) again using (22) it may 
be shown that (20) reduces to 


f+ic ^ 9(r) 


= - j- f b c f £ t ^ d 5 , ( b<r<c ) 


(25) 


/t -VF 


If we let 


t = p 2 , r = s 2 , /rg ( r ) = G(.s), /tf(t) = F(p) 


? 2 
b = d , c = e 


(26) 
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(25) becomes 


- 1 / , e d P 

7T d p “ s 


which can again be solved in closed form for a given G(s). [6]. 

For example, if the stamp has a rectangular profile, g(r)=0=G(s), 
and the solution (27) satisfying 


f^ c f(r)dr = - P 


(28) 


becomes 


f(r) = - i — , (b<r<c) (29) 

2tt [r ( /r- JE) ( /c-/r ) 

3. STRESSES AROUND THE APEX OF THE WEDGE . 

As indicated before, after obtaining f(t) by solving (20), 
the stresses and the displacement derivatives in the wedge may be 
evaluated by means of definite integrals having f(t), ( b < t < c ) as 
the density function. From the view point of fracture of the 
solid of particular interest is the cleavage stress a 06 (r,e) 
which is known to be singular at the wedge apex r=0 for wedge 
angles $ 0 >7T * The inversion of (11) would indicate that for small 
values of r the dominant term for the stresses is of the order 

j~r"^ s -l + 2 ^ , (i,j = r,9), 

D(s_ 1 ) = 0, Re (s_ 1 )<-l (30. a, b) 

where s -j is the first root to the left of Re(s)--1. By examining 
(30b) with D(s) as given by (9e) it may easily be shown that s -j 
i s real and 
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( 31 ) 


s n = - 2 for 0<9 <tt , 

- 1 o- 

-2<s_-j<-3/2 for ir<6 o < 2 tt. 

Thus, for 0 <tt the stresses at r=G are bounded and are of no par- 

o- 

ticular interest. For 0 >tt -they are singular at r=0, (30) gives 
the power of singularity, and the residue of the related inversion 
integral at s -| gives the asymptotic values. After routine manip- 
ulations from (11), (12), and (9) we find 

lime ^(gpe+IOgg) - - rwh • f b CtPf(t)dt ’ (32) 

r->0 “ • 

P = s _! +1 » 

D "( s_ 1 ) = 2[p(l-cos20 o ) - e o sin2p6 o 3» 

F(s-j) = pe“ l6p (c 1 + ic 2 )sih0 + ( c ] -i c 2 )e 1 9 si np6 , 

C 1 = p(l-cos26 o ) + (1-cos2pe o ), 

C 2 = psin2e o +sin2p0 o . (33) 

If we define 

2 + s_.j = a) , F ( s _ i ) = F -j ( 0 ) + i F 2 ( 0 ) * ( 3 ^ ) 

the "stress intensity factor" for the cleavage stress at r=0 may 
be expressed as 

F « ( 0 ) -i 

k(9) = lim r w a Q0 (r,0) = - — y / b t f(t)dt (35) 

r^O -1 

Thus, for a given wedge angle 0 Q the integral in (35) will be the 
measure of the intensity of stresses at the apex of the wedge. 
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4. THE CASE OF b=0. 


In the problem of flat-based stamp shown in Figure la and 
in the crack problem shown in Figure lb, the limiting case of 
b=Q is of special interest. In this case, as the crack problem 
would indicate, the singularity at the end point r=b=0 is that 
of a symmetrically loaded wedge of angle 20 q . In fact, it may 
also be shown that in the integral equation (20) part of the 
kernel corresponding to the inner integral is not bounded for 
all values of r and t in the closed interval [0,c], and becomes 
unbounded when r and t approach the end point r=0 simultaneously. 
This means that for b^O the kernel of the singular integral equa- 
tion (20) is of generalized Cauchy type and the related fundamen- 
tal function is of the following form: 

w(t) = t 6 (c-t) a s ( -1 <Re ( a , 0) <0) . (36) 

The character! s ti c equations giving the powers a and 3 are found 
to be 

cotTra = 0, 

(l + 3)sin20 o + sin2(l+3)8 o = O. (37. a, b) 

From (37) it is seen that a and 3 are real, a=-0.5, 3=0 for 
O < 0 o <'rr/2, and -1<3<0 for 2Tr>0 o >Ti/2. For example, 3 = - 0 . 5 for e o =Tr , 
3 = -2/3 for 6 *3Tr/2,.and 3^*1 for In the crack problem 

(i.e., 0 2 tt ) 3 = -l and it appears that the stresses have a nonin- 
tegrable singularity. Since the load is applied at the crack tip 
(i.e., b=0) this result is expected. 
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5. SOLUTION AND NUMERICAL RESULTS. 


The integral equation (20) subject to the condition (4) can 
easily be solved numerically once the stamp profile g(r) is speci- 
fied (for the numerical technique see, for example, [7,8]). Most 
of the results given in this section refer to flat stamp and crack 
problems shown in Figures la and lb, which are considered to be of 
greater practical interest. However, some results on the half- 
cylinder stamp problems shown in the insert of Figures 4a will 
also be given. 

In the crack and flat stamp problems, the problem may be posed 
in one of two ways. In the first one may specify the external load 
P and the distance of its line of application d. In this case, 
generally there would be a small rotation £ q of the stamp (or the 
half wedge away from the net section be). This quantity may be 
computed by using the following moment equilibrium condition: 

/ b c f(r)rdr = + Pd (38) 

where P is the magnitude of the applied load and noting that 
’f(r)=a £1Ci (r,0) , + sign refers to the crack and - refers to the 
stamp problem. 

In the second group of problems one may specify the load P 
and the angle of rotation e . In this case d is unknown and is 
again determined from (38). The practical problem here of course 
is the determination of d for a given load and no rotation 
(i.e., e =0). In both of these problems the input function in (20) 
is (see ( 3b) ) 
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(39) 


g (r) = |p u 0 (r,O) = e 0 , (b<r<c). 

In the flat stamp and crack problems for b>0 the fundamental 
function of the integral equation and the solution may be expressed 
as 

w(t) = [(c-tKt-b)]'* 5 , f(t) = F(t)w(t), ( b < t < c ) (40. a, b) 


where F(t) is bounded in b<t<c. In this case* in addition to the 
contact pressure - f ( r ) in stamp and the net section stress f(r) 
in crack problems, the quantities of physical interest are d or 
e , the stress intensity factors defined by 

k ( b ) = 1 i m /2 ( r-"b ) f ( r ) , k ( c ) = 1 i m /2 ( c-r ) f ( r ) , (41) 

r-*b r+c 

and the power to. and the strength k(0) of the stress singularity 
at the wedge apex for 0 >tt in stamp problems. 


Recalling the definition of stress intensity factor at r=0 

for e o > 7 r, (see (32-35)), we define a normalized stress intensity 

factor k by 
o J 


k(e) = lim r w (o 0 +i a QQ ) 


r+0 


re 


ee 


F 1 (0) + iF«(O) -j 

1 »-(s.p - V 1 "W 

F 1 ( 8 ) + i F 2 ( 6 ) c-b “- 1 

77 D ' ( s _ y ) 2 o' 


(42) 


In the case of b=0 the fundamental function of the integral 
equation (2,0) is given by (36) where a and 3 are determined from 
(37). In this problem too the solution may be expressed by (40b) 
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with F(t) being again a bounded function in 0<r<c. In this case 
the stress intensity factor at r= 0 is defined by 

'k(b) = k(0) = lim /7r' e f(r). (43) 

r^O 

Some calculated results for b>0 and e^-0 for geometries shown 
in Figures la and lb are given in Table 1. In this table the stress 
intensity factors k(c) and k ( b ) defined by (41) and (43) are pos- 
itive in the crack problem and are negative in the stamp problem. 

The normalizing factor for k(c) and k(b) are defined by 

k = — — , for b>0, 

n tt/I 

k n . = — y+q * for b = 0 , a = (c-b)/2 . (44. a, b) 

Tr.a 

The table also shows the powers of singularity m and 3 (see (43), 
(40b) and (36)). The values of w refer only to b>0 case where 

w=0 for O<0 <tt , 03=0.5 for 0 = 2ir, and 0<w<0.5 for ir<0 <2 tt. The 

o- o - - - o - 

values of 3 refer only to b=0 case. Note that as required by the 
physics of the problem 3{ 0 Q ) = - a> ( 20 q ) . Thus, the stress intensity 

measure k is nonzero only for co>0 (i.e., for 0 >7T,b>0), and k(b) 

o o 

is nonzero only for 3<0(i.e., for 8 o >tt/ 2, b = 0). 

For e Q =0 the values given for the distance d indicates that, 
as expected, d.+ (c + b)/2 as (c+b)+«> for a constant a=(c-b)/2, 
d>(c+b)/2 for 0 o <ir, d=(c + b)/2 for 0 o = tt and d<(c+b)/2 for 0 Q > tt , 
where d is the distance at which the resultant load P should be 
applied for zero rotation of the structure (Figure la) or the half- 
wedge (Figure lb). Another physically expected important phenom- 
enon may be observed from the results given for 0 o =6O°. It is 
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Table 1. The results for flat stam 

b> 0 s e =0 . 

-O' 


0 q ,0) ,8 

b/a 

c/a 

d/a 


10 

12 

11.277 

0 o =6O° 

8 

10 

9.335 

w=0 

6 

8 

7.424 


5 

7 

6.486 


4 

6 

5.5 69 


4 

6 

5.117 

0.=9O O 

2 

4 

3.192 

u * 
03=0 

1 

3 

2.279 


0.5 

2.5 

1 .858 

8=0 

0 

2 

1 .476 


2 

4 

3.037 

0 =120° 

1 

3 

2.056 

o 

II 

03 

0.5 ■ 

2.5 

1 .577 


0.2 

2.2 

1 .301 

8=-0. 38427 

0 

2 

1.136 

e n =i 50° 

0.5 

2.5 

mm 

S=o 

0.2 

2.2 


8=-0. 48778 

0 

2 

m 

V Plf 

QH 1 

0 

2 

1.0 


0.5 

2.5 

1 .496 


0.2 

2.2 

1.194 

BKlHiS 

0 ' 

2 

0.988 

mm 

1 

3 

1 .982 


0.5 

2.5 

1 .475 

■■■■■ 

0.2 

2.2 

1.165 

WBamnm 

0 

2 

0.939 


2 

4 

2.974 

6=270° 

1 

3 

1 .960 

30=0 . 45552 

0.5 

2,5 

1 .443 

0.2 

2.2 

1.120 

8=-2/3 

0 

2 

0.860 
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and crack problems, for 


k(c)/k n k { b ) / k n -k Q 


1.5486 0.4414 
1.6632 0.3225 
1.8356 0.1428 
1.9579 0.0151 
2.1183 -.1518 

1.2265 0.7583 
1.3634 0.5958 
1.5151 0.4003 
1 .6433 0.2238 
1.8436 ->0 

1.0684 0.9206 
1.1004 0.8731 
1.1319 0.8154 
1.1636 0.7350 
1.2005 0.8993 

1.0121 0.9817 
1.0153 0.9720 
1.0198 0.9958 

1.0 1.0 


0.9930 

0.9911 

0.9852 

1,0118 
1.0190 ■ 
0.4892 

2.821 

4.774 

0.9696 

1 .0446 

2.215 

0.9595 

1 1.0705 

2.862 

0.9483 

1.1166 

3.835 

0.9255 

0.4940 


0.9540 

1 .0592 

1 .779 

0.9323 

1.1007 

2.315 

0.9099 

1.1607 

2.917 

0.8842 

1.2716 

3.815 

0/8313 

0.7361 


































seen that at a value b=Ka, 4<K<5 the sign of the stress intensity 
factor k(b) changes (meaning that in the stamp problem the wedge 
"bends" and the contact stress around r = b becomes "tension" and 
in the crack problem it becomes compression). The technique de- 
scribed in this paper is also applicable for this case. The solu- 
tion of the crack problem is quite straightforward and may be ob- 
tained by using the formulation given in this paper provided one 
lets b=0 (i.e., the crack 0<r<b is now closed) and 3=0 at r=0, 
with the support of the integral equation (20) being (0,c). The 
stamp problem on the other hand, is somewhat different and has 
to be treated as a "receding contact problem" with the contact 

area c>r>b as an additional unknown determined from the condition 

o 

that the contact stress at r=b , (b >b) vanishes (see the last ex- 

o o 

ample given in this paper). 

Some sample results for the contact stress (or the net section 
stress) f(r) are given in Figures 2 and 3. In both cases the angle 

of rotation e is zero and the corresponding distance d of the 

o 

line of load application may be found in Table 1. Figure 2 shows 
the contact stress for 6 o =90°, b>0, and Figure 3 show the same 
for b=Q and for various values of the wedge angle 0 . 

The results for a symmetrically loaded stamp (i.e., d = (c+.h)/2) 
in the more interesting case of b=0 are shown in Table 2. In this 
problem one of the primary unknowns is the rotation e Q which is 
given in the third column. Here the normalizing angle e n is defined 
by 
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8 


P 



b/a r/a c/a 

Figure 2 Contact (or net section) stress for a 90-degree 
wedge (b^O) . 



Figure 3 Contact (or net section) stress for an elastic 
wedge (b*0) . 


Table 2. The results for the symmetrically loaded stamp or 
crack problem, b=0, d=c/2=a. 


o 

CD 

■ -3 

e c/ e n 

k(c) 

kn 

, °- 5+ e k(0 ) 

' aJ kn 

90° 

0 

-0.4727 

0.6486 


120 

0.38427 

-0.0967 

0.9132 

1.284 

240° 

0.66061 

0.0373 

1.0445 

0.449 

270° 

0.66667 

0.0829 

1.1010 

0.613 


te = (]+<) ? . (45) 

n 4pa 

The normalization factor k n for the stress intensity factors 
k ( 0) and k(c) shown in the table is given, by (44a) i.e., k n = P/TT/a. 
The results given in the table corresponds to the stamp problem 
where P is the magnitude of the applied (compressive) load, and 

e >0 means that the stamp rotates in the positive 0 direction 

o 

shown in Figure la. Only the results given for 0 o <it are valid 
for the crack problems. It should be noted that here even though 
one is dealing with the contact problem for a homogeneous medium, 

for b = 0 and- 0 >tt one obtains a stress singularity with a power 

o 

greater than 0.5 which is the strongest attainable power in notch 
and crack problems. 

The results for a semi-circular rigid stamp are given in 

Tables 3 and Figure 4. The normalizing stress intensity factor 

k which appearsin the table is defined by 
n • 

k n = /2 ( c'-FJ" 4 u/(1+k). (40 
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The constant k which is the measure of the strength of stress 

o 

singularity at the wedge apex for 0 is again defined by (42). 
In the example it is assumed that b=R=constant and the unknowns 
c and d are determined from the force and moment equilibrium, 
conditions (4) and (38). For the special case of 6 o ~tt the solu- 
tion is given by 


Table 3. The results for the semicircular rigid stamp, b=R. 


— “ — 1 '"I 

9 o 




■ 

n 






1.0338 

0.01115 

0.04596 


e =90° 


1 |f£ 

1 .0862 

0.06505 

0.1027 


0 


m§ 

1.1777 

■ 0.2379 

0.1752 






0.01187 

0.05055 

0.01166 

6 =270 



1 .0830 

0.07485 

0.1282 

0.07179 

0 


D 


0.3033 

0.2620 

0.2806 

0 = 1 80 

0 



1 .0333 

0.01178 

0.05 



P(r) 


- c-b 
" 1+K R 


(1 + 


>(££) 


r-b 
c-b' ' r-b 


(47) 


The resultant load P and its distance from the apex d are shown 
in Figures 4a and 4b, respectively. The results do not appear to 
be significantly different from those obtained for the half plane. 
Hence no extensive numerical work was carried out regarding this 


ex amp! e . 


The receding contact problem for small wedge angles is con- 
sidered as a last example. The results for a 60-degree wedge acted 
upon by a flat-ended rigid stamp are given in Table 4. The prob- 
lem has some anomalous features. In this case the index of the 
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singular integral equation is zero and the solution is of the 
following form 

f(r) = F ( r ) [ ( r-b )/ ( c-r ) , (b<r<c) (48) 


where F(r) is again bounded in b<r<c. In the numerical example, 
in order to avoid the iteration to determine b the problem is 
solved for. a constant stamp rotation, i.e.. 

Table 4. The results for a flat stamp on a 60-degree wedge, 
f(r)=F(r)I(r-b)/(c-r)]^. 


b/c 

a/c 

d/c 


mm 

mm 


0.7 

0.7025 

0.7050 

0.7075 

0.7100 

0.7104 

0.15000 
0.1 4875 
0.14750 
0.14625 
0.1 4500 
0.14480 

0.92401 

0.92466 

0.92530 

0.92595 

0.92660 

0.92670 


1 .619 
1 .627 
1 .634 
1 .641 
1 .649 
1 .650 

48.956 
64.708 
96.088 
189.27 
1 2339.8 

->-co 


g (r)=e o = constant , and. for each given b the load P is determined. 
The results are shown in Table 4. The normalization constants 
which appear in the table a re. defined by 


1 +K 

yc 


- y° 

l+K 


(49) 


A close examination of the results given, in the table would indi- 
cate that from the unloaded state P=0, e o = 0, b/c=0 (c=constant) as 
the load P is increased for a constant small rotation e Q b/c in- 
creases. It appears that there is a limit b/c=0.7104 beyond which 
the contact area cannot be reduced any further. At this value 
of b,e o -*0 for any fixed load P. From (20) and (28) with g(r) = E Q = 0 
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it is seen that the unknown function is f(r)/P and b=0. 7104c is 
really independent of the magnitude of the applied load P. This 
conclusion is typical of the general receding contact problems. 
From other trial solutions it was also found that it is not pos- 
sible to obtain any solution in the form of (48) (i.e., with a 
sharp stamp corner at c and smooth contact at b) for e Q <0. This, 
of course, is the result one would expect physically. 
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